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Abstract 

In this work we study the indirect effects of squarks and gkiinos 
via supersymmetric QCD radiative corrections in the decays of Higgs 
particles within the Minimal Supersymmetric Standard Model. We 
consider a heavy supersymmetric spectrum and focus on the main 
nondecoupling effects. We propose a set of observables that are sen- 
sitive to these corrections and that will be accessible at the CERN 
Large Hadron Collider and Fermilab Tevatron. These observables are 
the ratios of Higgs boson branching ratios into quarks divided by the 
corresponding Higgs boson branching ratios into leptons, and both 
theoretical and experimental uncertainties are expected to be mini- 
mized. We show that these nondecoupling corrections are sizable for 
all the proposed observables in the large tan /3 region and are highly 
correlated. A global analysis of all these observables will allow the ex- 
periments to reach the highest sensitivity to indirect supersymmetric 
signals. 
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1 Introduction 



One of the main goals of the next generation of coUiders will be to explore 
the Higgs sector phenomenology. The simplest candidates for Higgs sector 
physics beyond the standard model (SM) are the two Higgs doublet models 
(2HDM) and, among these, the leading one is provided by the minimal 
super symmetric standard model (MSSM) 0. This requires a Higgs sector 
called of type II, where one of the Higgs doublets couples to the toplike 
quarks and the other to the bottomlike quarks. The basic differences between 
the MSSM and a general 2HDM of type II (2HDMII) depend, first, on the 
supersymmetric (SUSY) sector, which is obviously absent in the 2HDMII, 
and, second, on the different values of the Higgs boson self-couplings. In the 
MSSM, because of the underlying supersymmetry, these couplings are fixed 
in terms of the electroweak gauge couplings. 

In all these models, the physical Higgs boson spectrum consists of two 
neutral CP-even Higgs bosons h° and H°, one CP-odd Higgs boson A°, and 
two charged Higgs bosons H"^ The tree-level parameters are the Higgs 
boson masses m/jo, muoi fn^o, and mH±, the mixing angle in the CP-even 
neutral sector, a, and the ratio of the vacuum expectation values of the two 
Higgs doublets, tan/3 = V2/V1. In the MSSM, due to supersymmetry, these 
can be written in terms of just two free parameters, which are usually chosen 
as mA° and tan/5. The Higgs boson masses are therefore not independent 
parameters in the MSSM, as they are in a general 2HDMII. Once the radia- 
tive corrections are included these tree-level mass relations are significantly 
modified, but there still remains a clear pattern for the Higgs boson masses 
in terms of the MSSM parameters Hence, if all these masses could be 

measured with good precision the mass pattern itself would be a first indirect 
indication of their SUSY origin. 

The most striking prediction of the MSSM is that the lightest Higgs 
boson mass lies below 130-135 GeV, within the reach of present and next 
generation colliders, the Fermilab Tevatron |^ and the CERN Large Hadron 
Collider (LHC) 0. The properties of this Higgs boson h°, however, are very 
similar to those of the SM Higgs particle in the so-called decoupling regime 
where mA° ^ 0- ^^e relevant h° production processes and decay 
channels at the LHC, this decoupling already occurs, in practice, at mA" 
values not far from the electroweak scale m-Ew- For instance, the difference 
in the width of the main decay h° bb, from the SM value is less than 
10% for mA° > 350 GeV and 2 < tan (3 < 50. There is just one region, 
corresponding to low m^o and large tan/? values, where the pattern of h° 
branching ratios significantly differs from the SM one. It is because in this 
region the dominant decays to bb pairs and to r+r^ pairs are enhanced and, as 
a consequence, the subdominant channels are significantly suppressed. Thus, 
for most of the (m^o, tan/?) parameter space it seems difficult to disentangle 
the nonstandard nature of h° and, in order to decide on its SUSY origin, one 
must consider in addition the phenomenology of the other Higgs particles 
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(for a review see, for instance, |^). The relevant question then becomes 
one of finding SUSY signals by looking into h° , H°, A° and production 
and decays at colliders, and of distinguishing these signals from those of a 
general 2HDMII. Obviously, if the SUSY particles are within the reach of the 
LHC and/or the Tevatron, the priority is to look for their direct production. 
However, it may well happen that these turn out to be too heavy to be 
produced directly. In that case one will have to look for indirect SUSY signals 
via their contributions as virtual intermediate states in standard processes, 
in much the same way as was done in the past at the CERN e~^e~ collider 
LEP for indirect effects from top and Higgs particles in precision observables. 
We are particularly interested here in their contributions to the radiative 
corrections involved in Higgs physics. 

In this work, we consider a scenario where some or all of the Higgs par- 
ticles have been discovered and their masses and branching ratios have been 
measured, but the SUSY particles have not shown up yet. We will study opti- 
mal strategies to look for indirect SUSY signals via their radiative corrections 
in Higgs boson decays [PHT5[|. Since it is known that there are specific decay 
processes where some SUSY radiative corrections do not decouple, even in 
the case of a extremely heavy SUSY spectrum |[T6|-p3|, we will concentrate 



on these particular channels. This nondecoupling behavior is a genuine part 
of Higgs sector physics and has not been found yet in other MSSM sectors, 
such as, for instance, electroweak gauge boson physics [^. The nondecou- 
pling processes we consider here are the Higgs boson decay channels into 
quark pairs, and more specifically we will concentrate on those that are en- 
hanced at large tan /5 values: h° bb, H" bb, A" bb ||,0,|TB]-[T|,gg,0 
llHTSl, . We will also analyze the top quark de- 

T^ , |l^ which is complementary to the tb 



and H+ -> tb 
cay channel t 



H+b 

decay for low m^o values. We will focus, in particular, on the SUSY QCD 
corrections from third-generation squarks, q, and gluinos, g, which are the 
dominant SUSY radiative corrections for most of the MSSM parameter space 
and are numerically sizable []9|-|T^,p!4|. For instance, in tb they can be 

as large as 60% and of either sign for quasidegenerate gluinos and squarks 
of mass 1 TeV and tan/5 = 50. These nondecoupling corrections were de- 
rived recently in for all the Higgs boson decay channels into all possible 
quark pairs, and analyzed in more detail for the h° — > bb [16^1^,2^,211 and 
— s> tb ||TP|, pO, p2[ channels. For nearly degenerate gluinos and squarks 



these SUSY QCD corrections, to one-loop level, can be generically written 



as r = r"(l + a5 



flMg 

Ml 



K)^ where F is the corresponding partial decay width 



without the SUSY QCD contribution, as is the strong coupling constant, /i 
is the bilinear Higgs boson MSSM parameter, Mg is the gluino mass, Msusy 
is a common generic SUSY mass characterizing the squark masses, and K is 
a function of (3 and a that depends on the specific channel p3[. The nonde- 
coupling behavior is seen here as a nonvanishing contribution to the partial 
width of order 0{asK), which is present even in the limit of a heavy SUSY 



spectrum where |/i| 



Msusy ^ ttiew- These corrections are known 
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to be particularly relevant for the h° , H°, A° —>■ bb, tb, and t H^b 

decays, because it is precisely in these channels that the function K grows 
linearly with tan/5 and can provide sizable contributions for the interest- 
ing large tan f3 region. It is worth mentioning that the SUSY electroweak 
(SUSY EW) corrections may also be relevant in particular regions of the 
MSSM parameter space p,[13|,0. In the large tan/5 regime and for quaside- 
generate squarks, neutralinos, and charginos, these are known to behave as 
~ {t^Y tan/?, where Msusy is the common generic mass for squarks, 

47r JKfsuSY 

neutralinos, and charginos, ht is the tree-level top quark Yukawa coupling, 

gmt 



s/2m.w sin/3 ' ^* ^'^^ quark trilinear coupling |[T6|-|T8|. For small 
or moderate At values and a heavy SUSY spectrum, these SUSY EW correc- 
tions are considerably smaller than the SUSY QCD ones and can be ignored, 
but for very large At values, \At\ ~ ~ Msusy ^ "^ew? these SUSY EW 
corrections do not decouple either and, depending on their sign, can increase 
or decrease the nondecoupling effect of the SUSY QCD corrections. 

In order to select optimal observables to look for these indirect SUSY 
QCD signals we require of them, first, to be measurable at the Tevatron [^], 
LHC 0, and/or the next generation linear colliders;^, second, to be most 
sensitive to the mentioned SUSY QCD corrections; and, third, to have the 
minimal theoretical and experimental uncertainties. We propose and analyze 
here a set of observables that satisfy all these conditions. These are the ra- 
tios of the branching ratios of the Higgs boson decays into third-generation 
quarks to the branching ratios of this same Higgs boson into third-generation 
leptonsQ We will show that these ratios can be considered optimal observ- 
ables because of the following reasons: 

• The SUSY QCD nondecoupling corrections contribute just to the numer- 
ator, that is, to quark decays and not to lepton decays. The first will be 
considered as the search channel, the second as the control channel. 

• These corrections are maximized at large tan /3 and are sizable enough to 
be measurable. 

• The production uncertainties are minimized in ratios. 

• They will be experimentally accessible at the LHC or Tevatron. 

• They will allow one to distinguish the MSSM Higgs sector from a general 
2HDMII. 

The paper is organized as follows. Section 2 is devoted to an experi- 
mental overview of Higgs particle searches at the Tevatron and LHC. The 
relevant regions in the (mAo, tan /3) plane are briefly reviewed. In Section 3 
the set of optimal observables is presented and analyzed in full detail. The 
leading contributions from tree-level and standard QCD corrections to these 
observables are also studied and their theoretical uncertainties estimated. 
The explicit formulas for nondecoupling SUSY QCD corrections from heavy 

"'^The case of linear colliders is not explicitly discussed here, but most of our results 
apply to them as well. This case has been considered recently in [ p5| . 
^Some preliminary results were presented by one of us in |26[|. 
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squarks and gluinos are presented. A discussion of the relevance of the SUSY 
EW corrections is also contained in Section 3. A comparison between a gen- 
eral 2HDMII and the MSSM is also included. The numerical results for the 
proposed observables and a discussion as a function of the MSSM parame- 
ters are presented in Section 4. Finally, the conclusions are summarized in 
Section 5. 



2 Experimental overview 



In this section we turn to the perspectives of the experimental measurement 
itself. We focus on the next hadron colliders: the Tevatron and LHC From 
this point of view, a measurement of the SUSY radiative corrections in Higgs 
physics is the next step in the problem of discovering a signal of the existence 
of one or several Higgs particles. Only after a discovery, but immediately 
after, will the question of what is the underlying structure of the Higgs sector 
be raised. As we have said, this is the principal motivation of our study. 

The discovery reach problem has been studied carefully and extensively 
by several working groups, both for the Tevatron Run 2 ^ and for the 
LHC 0. In particular, the searches for SM or MSSM Higgs bosons were used 
as benchmark channels during the design of the ATLAS and CMS detectors. 

Following the results of these studies, one learns that for values of tan /3 
large enough so that the corrections are sizable, there are three regions in 
the {mA°, tan /?) plane of phenomenological interest. 

The first is the region of large values of tan f3 (say > 10), and low m^o 
(say tua < 120 GeV). Here several production processes grow with tan^ /3, 
and the corresponding cross sections are much larger than the typical SM 
ones. One example is the associated production of a CP-odd boson A° and 
bottom quarks: gg{qq) — ^ bbA°. In this region of the MSSM, the lighter 
CP-even boson h° has essentially the same couplings and mass (within 5%) 
as the A", effectively doubling the production cross section [l|]. The heavier 
CP-even boson H° has SM bosonlike couplings, and its production rates are 
considerably smaller. In turn, the charged Higgs boson is lighter than 
the top quark, and t — > bH~^ happens in a sizable part of the top quark 
decays. 

With such large cross sections and relatively small masses this has been 
the main region examined at the Tevatron during Run 1. In fact, the Collider 
Detector at Fermilab (CDF) Collaboration has already searched for associ- 
ated production of 6M°(/i°), where A"{h") hi 0, or A"{h") t+t' 
excluding a part of the {mA°, tan/3) plane. Also, the CDF and DO Collab- 
orations have excluded a part of the MSSM parameter space corresponding 



to unobserved decays t — > bH~^, with ^ t^v [^, more specifically, the 
region where p^^^^+fc^^ ^1. It is clear that the search will continue in Run 
2, and a discovery is possible before the start of the LHC 0. 

The first analysis takes advantage of both the large cross section (0(1-100 
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pb)) and large 3(^4° 66) ~ 90%. The experimental signature is four jets 
with at least three of them b tagged. This channel has also the advantage 
that it allows the reconstruction of m^o and, once this mass is known, the 
production rate gives an idea of the size of tan f3. On the other hand, as 
in any search in a purely hadronic channel, the QCD background has to 
be properly understood. Systematic errors at the 20 — 30% level are not 
infrequent ^ . 

In contrast, the channel A" — » r+r^, has to overcome a relatively smaller 
branching ratio (~10%). However, the electroweak background > r+r~) + 
jets is considerably less dangerous . 



If an A° signal were found in these analyses, the corresponding ratio of 
bb to r+r^ events would be a direct test of the SUSY QCD corrections. 

A possible signal in the charged Higgs boson channel could also be use- 
ful |]T3[. In this case, it also happens that the SUSY QCD corrections do 
not decouple in the t — > bH'^ decay but they do decouple in the standard 
channel, t bW^ |^,^. This is the reason why r(t bW'^) is considered 
now as the corresponding control width. The problem in this case is that 
the tree-level prediction changes rapidly with mH+ and tan (3. On the other 
hand, the accuracy in the experimental reconstruction of mH+ in the decay 
— > r+z/ is limited, and this channel by itself does not provide information 
about tan (3. It should be pointed out that the LHC experiments will be able 
to explore this channel very efficiently. 

The second region of interest still has large tan /3, but 140 GeV < ruA" ^ 
500 GeV0. With such large masses this can be considered the genuine search 
zone of the LHC experiments [^]. In this zone, the heavier CP-even boson H° 
has similar couplings and mass than the A° boson, again doubling the cross 
section of the process gg{qq) bbA"{H") (C(10-1000 pb)). Now the light 
CP-even Higgs boson h" behaves like the one in the SM with rriho ~ 110 — 130 
GeV (the precise mass value depends on the choice of the MSSM parameters). 
In this region, the mass of the charged Higgs boson is either very close to 
rrit and B(t — > bH^) is negligible, or mH+ > and t bH^ is directly 
forbidden. Then, the main possibility of producing charged Higgs bosons is 
in the process gg{qq) — »■ btH^ . Again, this cross section grows approximately 
(at large enough tan/?) with tan^ (3. 

The main possibility of the Tevatron here would be to detect the SM-like 
boson h in the associated production qq W{Z) h° and h° — * bb. Also, the 
analysis of the ratio {A° bb)/{A° t^t^) can be extended, but only for 
relatively small values of mA° and very large values of tan {3 0. 

With respect to these last channels at the LHC, the challenge will be to 
detect the decay A°{H°) bb in events with four b jets. Studies have shown 
that to control the QCD background will be very difficult. In contrast, it will 
be possible to measure the decay A°{H°) t^t~ and even A°{H°) —>■ fi^fi^ . 

^The properties of the Higgs sector in the intermediate region 120 GeV < mA" 140 
GeV are a mixture of the two cases discussed. A detailed description is complicated and 
beyond the scope of this paper. 
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Both channels are leptonic, and therefore do not allow an estimation of the 
SUSY QCD corrections. However, they will provide good measurements of 
m^o and tan /3, which can be used in addition to other analyses (say, — > 
T^v). Current studies for the LHC indicate that values of Am^o/m^o 
better than ~ 1% and of A tan /5/ tan/? better than ~ 12% ( if tan/3 > 10) 
can be reached. 

In the case that m//+ > mt + mt, the decay channel th opens 

and can be used not only to complement the decay — > r+z/, but also to 
provide an independent estimation of mA° (in the frame of the MSSM) and 
tan /?. 

Finally, we mention the very interesting region with mA° ^ 140 GeV and 
intermediate 3 < tan /? < 10. At the LHC this zone is called the "hole" 
because the only Higgs particle accessible is the SM-like one. Here the main 
question would be not only to discriminate the MSSM from a more general 
2HDM sector, but to know if there is something beyond the plain SM. 

We will see in the next sections how the SUSY QCD corrections can be 
used for that matter. 

3 Optimal observables 

In this section we present the set of observables that are the most sensitive 
to the nondecoupling SUSY QCD corrections. 

Since we are looking for indirect signals of SUSY in an experimental sce- 
nario described by the discovery of Higgs particles, with masses following the 
MSSM pattern, and where the SUSY particles are too heavy to be produced 
in colliders, we will search for observables with very specific conditions. From 
the theoretical point of view the requirements are the following: 

• In order to discriminate between the MSSM and a nonsupersymmetric 
2HDMII, we need observables with different predictions for these two 
models. 

• These predictions should be distinguishable even in the case of a very 
heavy SUSY spectrum. 

• The theoretical uncertainties must be minimized. Some specific ratios 
of branching ratios will cancel these uncertainties either totally or par- 
tially. In particular, we wish the tan (3 and a dependence to appear 
just in the correction, but not in the lowest order contribution. 

On the other hand, the experimental requirements are the following: 

• A control channel is needed in order to eliminate systematic errors, so 
that ratios of event rates are better than event rates themselves. 

• The corrections to the observables should be sizable in the large tan (3 
regime, since the relevant Higgs particle production cross sections grow 
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with this parameter, such that we have better statistics for larger values 
of tan 13. 

• There should be good identification of the final state particles. 

• Experimental uncertainties from Higgs particle production must be 
minimized. These will cancel in some specific ratios of events from 
Higgs particle decays. 

• The expected accuracy at the LHC and Tevatron for these observables 
should be good enough that the deviation produced by the corrections 
becomes apparent. In particular, the size of the SUSY QCD corrections 
must be larger than the expected error bars. 

We propose the following set of optimal observables, which satisfy the 
previous theoretical and experimental requirements: 



Bih^^bb) _ B{H°^bb) 



B{h° T+T-) ' B{H° T+T- 

_ BjA^^bb) _ BiH+^tb) 

= -rn-Tn ' = 



B{A° T+T-) ' " B{H+ T+u) 



In addition, we consider the following ratio: 

^ B{t Hn) 
B{t^w+by ^ > 

which complements the charged Higgs boson observable in the low m^o re- 
gion. The predictions for these observables at the tree level are the same for 
any general 2HDMII, and in particular for the MSSM. These are 



Otree at ^bf^b i 1,0 tto . ,otree at '^bf^b 

_ 2 - - mg) (m^ cot^ [3 + ml tan^ (3) - Amjml ^ 1/2 

'^11+ — J-^rlllu+ T. r, „ , 2 ^^^^^^ ^ 



O 



tree _ (mg + - m\^){m\ cot^ /3 + tan^ ji) + 4mgm,^ 



(3) 



where 

= ^1 - 4m2/M2 , X^l = ^(1 - {ml + ml)ymt){l - {ml - mir/mt). 

As a general remark, the differences in the predictions from the various 
models for all these observables come in the corrections beyond the tree 
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level and will be discussed in the following subsections. In particular, it is 
interesting to emphasize that the predictions for the neutral Higgs boson 
observables at the tree level are independent of tan/3 and a. The charged 
Higgs boson observable depends on tan /3, but this dependence is very mild in 
the large tan /3 region. Therefore, these observables are especially sensitive to 
any extra contribution growing with tan /3. In the top quark decay observable, 
however, the tree level prediction goes as tan^ /5, for large tan /5, and it will 
be more difficult to disentangle the mentioned contributions. Other analyses 
of the observables for the neutral channels in the large tan (5 limit and in the 
zero external momentum approximation can be found in The neutral 

channel case with Higgs boson mass corrections included has been analyzed 



m 



3.1 Supersymmetric QCD contributions 

The nondecoupling SUSY QCD contributions to these observables in the 
MSSM can be easily derived from the results for the effective Yukawa inter- 
actions of the Higgs sector with top and bottom quarks. These have been 
recently obtained at one-loop level, for arbitrary tan /5, and by a functional 
integration of bottom and top squarks and gluinos in . For nearly degen- 
erate heavy squarks and gluinos, these corrections can be written as 



Oh° = O'^o 
Oa° = 0% 



Ot =01 



2as 




3tt 




2as 




3tt 


^^-'SUSY 


2as 




37r 


^*^SUSY 


2as 




Stt 


^^-'SUSY 


2as 




371 


^'^SUSY 



(tan [3 + cot a) 
(tan/3 — tana) 
(tan 13 + cot /3) 
(tan [3 + cot 0) 
(tan j3 + cot j3) 



(4) 



Here, /i and Mg are the MSSM bilinear parameter and the gluino mass, 
respectively. Msusy is the common SUSY mass for squarks and gluinos 
(MsusY ~ ~ Mg) and as is the strong coupling constant evaluated at 
the corresponding decaying particle mass. The leading terms 0° refer to 
the value of the observables without the SUSY particle contributions and 
will be discussed in subsection 3.3. These formulas are for heavy squarks 
and gluinos, that is, for large Msusy, and are valid for all tan/3 and mA° 
values. Corrections to the previous formulas are either of higher order in as 
or suppressed by higher inverse powers of the heavy SUSY masses Msusy 
and, for the present analysis, can be safely ignored. 
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From a first look at eq.(^ one sees clearly that the observables have 
different predictions within the MSSM and in a general 2HDMII, since the 
former includes the SUSY QCD corrections, which do not exist in the latter. 
This main difference comes from the nondecoupling behavior of the SUSY 
QCD corrections in the Higgs boson decays into quarks (and in the top quark 
decay into charged Higgs bosons), which is relevant even for a very heavy 
SUSY spectrum, such that Mg ~ ~ Msusy ^ ?^ew- Notice also that, as 
mentioned before, these SUSY QCD contributions grow linearly with tan/5 
so we expect sizable corrections for large tan/5. 

There are two limiting situations that are worth mentioning: the large 
mA° limit m^o ^ rriz, and the large tan/3 limit tan/5 ^ 1. In the former 
limit, cot a approaches — tan/5 (correspondingly, tana —>■ — cot /5) and, as 
can be seen in eq.(^), the SUSY QCD correction decouples in Oho. The other 
observables get exactly the same nondecoupling correction, proportional to 
(tan /5+cot /5). In addition, in this limit, the H° and masses approach the 
large A° mass. The leading contribution O^o approaches the SM prediction 
and, therefore, as noted before, the situation cannot be distinguished from 
the SM case. In the large tan (3 limit, the corrections do not decouple in any 
channel, they are all proportional to tan/5 and all have the same sign. The 
universal character of the corrections will be very useful in a global analysis, 
because it will yield strongly correlated signals in all the channels. 



3.2 Other nondecoupling contributions in the minimal 
supersymmetric standard model 

It is well known that, in particular regions of the MSSM parameter space 
the radiative corrections from the SUSY EW sector can also be relevant 
in the decays of the Higgs bosons and top quark that we are considering 
here [|], |13|, |1^] . These SUSY EW corrections are largely dominated by the 
contributions with charginos or neutralinos and third-generation sfermions 
in the loops; the contributions from the Higgs sector being negligible. For 
instance, in the tb case with 171^^+ = 250 GeV and for light top 

and bottom squarks and light charginos and neutralinos, Mg, My. < 200 
GeV, they contribute as much as ±20% with respect to the tree level width, 
whereas the SUSY QCD loops induced by squarks and gluinos are by far the 
leading SUSY effects and give a contribution larger than ±50% |13|. These 



SUSY EW corrections can only be competitive with the SUSY QCD ones in 
the particular region of the SUSY parameter space where tan /5 is very large 
(tan/5 > 20), the bottom squarks are very heavy {M^ > 300 GeV), and the 
top squarks and charginos are relatively light (M^, ~ 100 — 200 GeV). For 
this particular choice, the total SUSY correction remains around (30-50)% 
of the tree level width, with at most half of it of SUSY EW origin | |13[| . 

For our present assumption of a very heavy SUSY spectrum with nearly 
degenerate SUSY particle masses (iW^ ~ Mj^ ~ Mg ~ M^ ~ ~ Mgugy ^ 
mEw) the SUSY QCD corrections provide by far the dominant contribution 
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to the one-loop MSSM correction. However, it is worth analyzing in more 
detail the SUSY EW corrections, which are known to provide extra non- 
decoupling contributions. These nondecoupling SUSY EW effects have been 
estimated in the literature just in the large tan /3 3> 1 regime and by means 
of an effective Lagrangian formalism, which works in the zero external mo- 



mentum approximation |jT6|, |T8|, [T9| . In this approach, the potentially large 
tan (3 enhanced SUSY corrections to the Higgs boson-quark-quark Yukawa 
couplings are induced via the quark mass corrections. In particular, the ex- 
pression for the bottom quark Yukawa coupling, at the one-loop level and in 
the large tan /? limit, is |T|,|19|,|g] 

hf, = — tan/5(l — Amfe), (5) 

V 

where rrih is the on-shell pole bottom quark mass, v = 174 GeV, and the 
tan (3 enhanced radiative corrections are encoded in 

Am, = Amf«^^ + Amf^^ (6) 

where Am^*^*"^ and Am^^^ refer to the bottom quark mass corrections 
from the SUSY QCD and SUSY EW sectors, respectively. Am^*^*"^ is dom- 
inated by the bottom squark-gluino loops, and to leading order in the strong 
coupling is @,|l6|,|l|,|g: 



^^SQCD ^ ^M^/xtan/3/(M^^,M^^,M5). (7) 

For sizable values of the soft trilinear coupling At, Am^^^ is dominated by 
the top squark-chargino loops and more precisely by the top squark-charged 
higgsino contribution. By neglecting the bino effects, which have been found 
to be numerically insignificant, the SUSY EW mass correction to leading 
order in the top quark Yukawa coupling and the electroweak gauge coupling 
is i,H,|l], 



1,2 2 

^"^^'"'^ = T7r^Afttanpi{Mi^,Mi^,ix) - J^M^^itanP [c2j(M^^, M2, /i) 

where M2 is the SUSY soft breaking chargino mass parameter and Cq = cos 6q, 
Sq = sin 9q, with 9q the q mixing angle. The one-loop integrals in Am^'^'"^ 
and Am^^^ are defined by, 

2 2 r ™? I 2 2 r 1^0 , 2 2 7 "i? 

mimiloq^ + mimiloq^ + mimiloq^ 
/(mi,m2,m3) = 3-—^ ^ . (9) 

For the simplest assumption that is considered in this work of all SUSY 
soft breaking parameters and the fi parameter being of comparable size. 
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<=?~^~i + oh^ . (10) 



Mq ~ Mfj ~ Mj^ ~ Mg ~ Mi^2 ~ ~ l/^l ~ MsusY ^ "^Bvy, the loop 
integrals and mixing angles behave as, 

I{m^,m2,m3) ~ — 2 + 0[ —3 , , 

' 2^^[Msusy 
and consequently the mass corrections behave as, 

A SQCD Ois MgfJ, 

First, we see here that the nondecoupling effects induced from the previous 
expression for Am^*^*"^ on the bottom quark Yukawa coupling through eq. (|^) 
are exactly the same as the ones extracted from the observables Oj;/ in eq.(^) 
in the large tan (3 limit. Therefore both approaches coincide in this limit, but 
just eq.(|[) gives the correct behavior for moderate and low tan/3. Second, we 
see that, even in the case of very large tan/5 values and large /i. At and M2, 
the size of the SUSY EW corrections remains always well below the SUSY 
QCD corrections. The relative signs of these contributions depend on our 
choice of the relative signs of Mg, /i. At, and M2. By choosing a combination 
of signs that maximizes the size of Am'^^^ and for equally large SUSY 
parameters we find that lAmf'^^l is less than 50% of \Am^'^'"^\. Therefore, 
for this assumption on the size and signs of the SUSY parameters, we can 
conclude conservatively that the effect on the observables of eq.(^ would be 
reduced by at most 50%. 

Since we are studying here the sensitivity of the observables in eq.® 
involving Higgs bosons and top quark decays to the SUSY QCD nondecou- 
pling corrections for all tan j3 values, a more precise analysis of these SUSY 
EW nondecoupling effects would require the use of compact formulas valid 
for all tan/? values, which are not available so far in the literature, and it 
is beyond the scope of this work. We will, however, take them into account 
conservatively in the final numerical analysis and conclusion by considering 
this somewhat pessimistic case in which they conspire for a 50% reduction 
of the SUSY QCD signal. 



3.3 Predictions for 0° in the minimal supersymmetric 
standard model 

Here we discuss the predictions for the leading contributions to the observ- 
ables, 0°, that do not include the SUSY particle contributions. Our purpose 
is to evaluate their uncertainties, and study whether or not they can mask 
the SUSY QCD corrections. Uncertainties below 1% are neglected. 
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The leading contributions 0° are computed here within the MSSM, and 
consist of the tree level part plus the one-loop as corrections from standard 
QCD0. They are evaluated numerically with the FORTRAN program HDE- 
CAY [^], using the renormalization group approach of ||3^ to obtain the 
MSSM Higgs boson masses. These standard QCD corrections contribute ob- 
viously just to the decays into quarks and are known to be as large as 50% 
(for a review, see In order to take into account large contributions from 
higher orders, we have used the running instead of the pole bottom quark 
mass. This resums the leading logarithms and improves the convergence of 
the series. To illustrate this, we show here the case of A° — >■ bb. We compare, 
in Table 1, the results of r{A" bb) for m^o = 500 GeV (normalized to the 
tree level width), computed at tree level, and at orders 0{as) and C(a|) in 
perturbation theory, using both the pole and running masses. 





tree-level 


0{as) 


0{al) 


Pole mass 


1 


0.523 


0.369 


Running mass 


0.309 


0.351 


0.363 



Table 1: Comparison of QCD contributions to r{A" bb), normalized 
to its tree level value, using the pole bottom quark mass and the running 
bottom quark mass, for m^o = 500 GeV. 



As can be seen from the results in this table, the convergence of the series 
is notably improved when the running mass is used. In addition, the error 
committed by ignoring C(«|) corrections is reduced to values ~ 3%. 

In the numerical evaluation of all the 0° we will include, therefore, these 
one-loop 0{as), QCD contributions and use the running bottom quark mass. 
Both the running bottom quark mass and as are evaluated at the corre- 
sponding particle decaying mass. We do not include, however, the standard 
electroweak radiative corrections from W^, Z, 7, or extra Higgs bosons since 
they are known to be below ~ 1% in the MSSM (for a review, see p| for the 



Higgs bosons decay case and |jT5[ for the top quark decay case). 

The predictions for 0° as a function of mA° and tan (3 are shown in Fig. |l[ 
The contour lines for 0° in the {m a°, tan f3) plane show clearly a very mild 
dependence on tan f3, for large tan (3 values, in all the Higgs boson channels. 
In the neutral channels, this comes exclusively from the dependence of the 
MSSM Higgs masses on tan /3. The dependence on m^o comes mainly from 
the running of the bottom quark mass, which is evaluated at the correspond- 
ing Higgs boson mass. In the case of h°, this dependence is frozen in the large 
mA° region, where it behaves like the SM Higgs boson. In the case of the 
H~^, the prediction grows significantly with m^o, as mH+ separates from m^. 
Finally, for the top quark decay channel, the prediction changes rapidly with 



*This is to be consistent with the SUSY QCD as corrections. 
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mA° and tan j3. The tan^ (3 dependence can be seen in the large tan (3 region. 
Moving to the right part of the plot, the size of the predictions decreases as 
mH+ approaches m^. 

We next analyze the theoretical uncertainties in the evaluation of the 
ratios in eq. (|1[) and eq. coming from the experimental errors in the values 
of the SM parameters involved in their determination, at one-loop and order 
as- We have found that the only significant errors are the ones coming 
from errors in the bottom and top quark masses and in the strong coupling 
constant as- We have used the results from [p^] : 

= (4.6±0.2) GeV, rrit = (174.3±5.1) GeV, as(M|) = 0.118±0.002 

The uncertainty coming from these experimental errors is 60 = q'^ , 
where O is the observable evaluated taking central values for m^, rrit, and 
as, and O' is the observable evaluated with one parameter reaching one 
extreme value. In Table 2 the different uncertainties are summarized for 
three different values of mA°, in the region tan/? = 10 — 60. We use 60n = 
60 h", 60 H°, 60 A", and 60c = 60 1 for m^o = 100 GeV, while 60c = 60 h+ 
for TUAo = 250 and 500 GeV. 

Notice that the corrections are approximately constant as a function of 
tan i3- The uncertainty related to the error in rrib is the dominant one, and of 
the same size for all the observables. Finally, the uncertainty related to rrit 
is significant for the charged Higgs boson and top quark cases, only in the 
vicinity of the corresponding kinematical thresholds. The total uncertainty 
in the sixth column has been evaluated as the sum in quadrature of the 
numbers in the previous three columns and the error coming from neglecting 
the corrections. 



6 


mA" 


as 




rrit 


Total 




100 GeV 


2% 


11% 


< 1% 


12% 


60n 


250 GeV 


3% 


10% 


< 1% 


12% 




500 GeV 


3% 


11% 


< 1% 


12% 




100 GeV 


1% 


8% 


6% 


11% 


60c 


250 GeV 


1% 


10% 


7% 


13% 




500 GeV 


3% 


11% 


< 1% 


12% 



Table 2: Uncertainties in the ratios induced by the experimental errors 
in rrit, ^b, and as for three values of rriA", in the region tan/3 = 10 — 
60. We use 60n = 60h°,60H°,60A° and 60c = 60t for rriA" = 100 
GeV, and 60c = 60 h+ for mA° = 250, 500 GeV. The total uncertainty 
has been evaluated in quadrature and includes the error coming from 
neglecting the corrections. 

In summary, the MSSM predictions for the previous Higgs boson ob- 
servables are known up to an uncertainty of the order of 13%. Any extra 
correction must be larger than that to be visible. 
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Figure 1: O" contour lines in the (mA°, tsji P) plane predicted within the 
MSSM. The input parameters are defined in the text. 
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3.4 Comparing a type II two higgs doublet model and 
the minimal supersymmetric standard model 

Here we investigate the similarities and differences of the set of observables 
defined in Section 3.1, in a general 2HDMII as compared to the MSSM. For 
a given Higgs sector mass pattern that is compatible with the MSSM, the 
values of the observables at the tree level, and after including the standard 
QCD corrections, coincide with the corresponding ones in the MSSM. Note 
that the mixing angle in the neutral sector, a, is a derived quantity in the 
MSSM, whereas in the 2HDMII it is an independent parameter (in principle 
— 7r/2 < a < vr/2, but actually the range in a is further restricted if the Higgs 
potential is required to be bounded from below) and, hence, its value can 
be different in the two models. However, since the branching ratios in the 
numerator and denominator of our observables O have the same dependence 
on a, at the tree level the O*''*^'^ themselves are a independent (see eq.(0)). 
When the standard QCD corrections are included, this a independence is 
still maintained. On the other hand, the standard electroweak corrections 
from IV^, Z, 7 are also similar in both models and, as already said, they 
have been estimated to be below 1% and can be ignored here. Therefore, 
the only differences, apart from the SUSY corrections, could come from the 
Higgs sector radiative corrections. As mentioned before, these are negligible 
in the MSSM but, in principle, could be larger in the 2HDMII. The reason 
is that these corrections involve the Higgs boson self-couplings which in the 
2HDMII are not restricted by SUSY and, in principle, can be large. We have 
estimated these potentially different contributions from the Higgs sector to 
the one-loop level and for a Higgs boson mass pattern compatible with the 
MSSM, and they turn out to be very small, certainly well below the common 
ones discussed in the previous section. 

In order to illustrate this we discuss the case of A° . We show in Fig.^ 
the maximum absolute value of the difference between the MSSM and 2HD- 
MII Higgs sector corrections for 0\o (after scanning the allowed range in 
a), relative to the tree- level prediction O^o*^. More specifically, these differ- 
ences come exclusively from the a-dependent one-loop diagrams, and from 
the triangular one-loop diagrams involving the Higgs boson self-couplings. 
The contour lines separate the regions in the (m^o, tan /3) plane where the 
maximum difference between the MSSM and 2HDMII is larger than 0.1%, 
between 0.05% and 0.1%, between 0.01% and 0.05%, and below 0.01%. In 



the region of interest, mA° > 84.1 GeV [0, this difference is always less than 
0.1%. 

In summary, as announced, the basic difference between the MSSM and 
a general 2HDMII model, with respect to our observables, is that the large 
SUSY particle loop corrections do not exist in the second case. In the next 
section we discuss these contributions. 
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Figure 2: The maximum absolute value of the difference between the MSSM 
and 2HDMII Higgs sector corrections for 0\o (after scanning the allowed 
range in a), relative to the tree-level prediction, O^r. Note that the area 
uiA" < 84.1 GeV (if tan P > 1) is already experimentally excluded at 95% 
CL Ml. 



4 Results and discussion 

In this section we present the numerical results for the SUSY QCD corrections 
to the set of observables Oho, Oh", Oa", 0^+, and Oj. The relevant param- 
eters in this discussion are tan/5, m^o, and Msusy (^susy ~ Mq ~ Mg). 
We first discuss the behavior with tan/5. In order to focus on this tan/3 
dependence, we consider the simplest choice for the SUSY mass parameters, 
that is, Msusy = = \fi\. We show in Figs. the results for three 
mA° values in the large, medium, and low m^o regions: m^o = 500, 250, 
and 100 GeV, respectively. The central lines in these figures follow the pre- 
dictions for the observables without the SUSY QCD contribution, namely, 
0°. The corresponding total theoretical uncertainties, discussed in Section 
3.2, are shown as shaded bands around the central values. The bold lines 
represent the SUSY QCD corrected predictions for > (solid) and /i < 
(dashed). The observable Ot appears in the lower left plot of Fig. ^ replacing 
Oh+, since in this case rut > mH+ + rrif,. 

One can see from the figures that, first, the predictions for all the observ- 
ables separate from the central values as tan P grows. The sign of the SUSY 
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10 10 

tan /? tan (3 



Figure 3: Predictions for Oh", Oh°, 0^+, and Oa° as a function of tan (3, 
for m^o = 500 GeV. The central lines are the corresponding predictions 
for 0° . The shaded hands cover the theoretical uncertainties estimated in 
Sect. 3.3. The hold lines represent the SUSY QCD corrected predictions for 
MsusY = Mg = |/i|. 



QCD corrections is positive for /i < and negative for /i > 0. The central 
values and their theoretical uncertainties are rather insensitive to tan (3 in 
the H° and A° channels. For the h°, there is a very slight dependence at 
low tan j3, while for the this dependence is very strong in the low tan (3 
region and it softens for tan/? > 15, as expected from eq. (3). In any case, 
the predictions for all the Higgs boson observables without the SUSY QCD 
contributions are tan/3 insensitive in the large tan/? region. This is very 
different for the top quark observable: it depends strongly on tan (3 in all the 
parameter space. 

Concerning the size of the SUSY QCD corrections, we see that there is 
always a tan f3 region where these are larger than the theoretical uncertain- 
ties. For mA° = 500 GeV and m^o = 250 GeV, the observables Oa" and 
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Oh° behave similarly and the predictions are outside the shaded band for 
tan/3 > 5. For mA° = 100 GeV, the crossing still happens at tan/3 > 5 for 
A° , whereas for H° the bold lines lie outside the error band in the whole 
region 2 < tan/? < 50. 

In the h° case the situation is qualitatively different. For mA° = 500 GeV, 
the SUSY QCD correction is below the theoretical uncertainty for all values 
in the region 2 < tan/3 < 50. This is a manifestation of the decoupling 
of this correction for large m^o values, as commented in Section 3.1. For 
smaller values of my^o, such as m^o = 250, 100 GeV, the decoupling has not 
effectively operated yet and the SUSY QCD corrections are sizable for large 
enough tan/5. In particular, for m^o = 250 (100) GeV, they are larger than 
the theoretical error band for tan/5 > 15 (5). 

Regarding the charged Higgs boson case, two different situations must 
be considered. For m^o = 100 GeV, where the decay into a top and a 
bottom quark is not kinematically allowed, the observable Ot is considered. 
Otherwise, for mA° = 250 GeV and mA° = 500 GeV, the relevant observable 
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10 10 

tan /3 tan jS 



Figure 5: Same as in Fig. ||, but for uiao = 100 GeV. The prediction for 
Ot appears in the lower left plot instead of 0^+, since in this case rrit > 



is Oh+ ■ As can be seen in Fig. |], the SUSY QCD corrections in Ot for mA° = 
100 GeV are above the theoretical uncertainty for tan/5 > 15. However, as 
pointed out in Sect. 2, the central value prediction also depends strongly 
on tan (3 (and m^o) and it will be difficult to identify the effect of the SUSY 
QCD corrections above the additional uncertainty related to the experimental 
errors on the measurement of tan/? and m^o. For Oh+ and m^o = 500,250 
GeV both requirements, the correction being larger than the theoretical error 
and the leading contribution being insensitive to tan/3, are satisfied for values 
larger than about 15. 

We next discuss the Msusy dependence. For this purpose, we fix the 
value to 250 GeV and choose Mg = Msusy- Figure ^ shows the predictions for 
the observables as a function of Msusy for mA° = 250 GeV, tan/5 = 30, and 
both signs of the fi parameter. Since, in this case, the SUSY QCD corrections 
vary as I/Msusy? their size will be below the theoretical uncertainty for large 
enough Msusy- This occurs in Fig. ^ for Oh", Oa", Oh+, and Ot at about 
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1500 GeV. For Oh" it is at about 500 GeV. Therefore, except in the hght Higgs 
boson case, we expect the observables to be sensitive to indirect signals of 
super symmetry via these SUSY QCD corrections up to quite large values of 

^SUSY- 

We finally turn to the question of the experimental resolution required to 
evidence the above effects, after a Higgs boson signal is found at a given posi- 
tion in the (m^o, tan /3) parameter plane. The predictions for the observables 
including the SUSY QCD corrections for Msusy = Mg = are shown in 
Fig. 0for /i < and in Fig. ||for fi > 0. Comparing these figures with Fig. ^ 
the patterns of the O (solid) contour lines change noticeably with respect to 
the ones for 0°. The differences are mainly because of both the large size of 
the corrections and the different behavior with tan (3. 

The shaded zone in each plot represents the region where the corrections 
(if existing) are totally hidden inside the theoretical uncertainty discussed 
in Sect 3.3. Even the perfect experiment with infinite statistics and perfect 
resolution cannot make conclusions about the size of the corrections. This 
zone is particularly large for the SM-like h° covering a large fraction of the 
LHC "hole" (see Sect. 2). On the other hand, to have a decoupling channel 
can be useful as it can provide the "calibration" of what should be expected 
for the other channels without corrections. 

In all plots in Figs. |^ and ||, the space above the long dashed line is the 
plane region that could be accessed experimentally with a modest resolution 
of 50% in O. Except for the h°, nearly all points with tan/5 > 20 — 25 
are testable at the one sigma level. The zone above the short dashed line 
requires an experimental resolution of 20%. Again with the exception of h°, 
this zone covers approximately the region tan/3 ^ 10 — 15. Therefore, if an 
experimental resolution of 20% can be achieved at the Tevatron and LHC, 
the analysis of the observables proposed in Sect. 3 could be used to search 
for indirect signals of SUSY in the main part of the relevant regions 1 and 2 
discussed in Sect. 2. 

Finally, we estimate the effect of the extra non-decoupling corrections 
from the SUSY EW sector. As discussed in Sect. 3.2, for very large At val- 
ues, the SUSY EW corrections do not decouple and this could increase or 
decrease the nondecoupling effect of the SUSY QCD corrections. For the 
present assumption of equally large SUSY parameters and by choosing a 
particular combination of signs that maximizes the size of Am^^^ , we find 
that this contribution always remains below 50% of the Am^*^*"^ one. In 
order to be conservative, we have performed the exercise of considering this 
pessimistic case where the SUSY EW corrections reduce the SUSY QCD sig- 
nal by 50% (see Figs. ^ and [l^). As in Figs. ^ and ||, the shaded zone in each 
plot represents the region where the corrections are completely hidden inside 
the theoretical uncertainty. The long (short) dashed lines again limit the 
zones where an experimental resolution of 50% (20%) is required to achieve 
a meaningful measurement. 

Figures ^ and ^ show that, as expected, the reachable values of tan/3 
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Figure 6: Predictions for Oh", Oh°, Oh+, Oa", and Ot as a function of 
MsusY = Mg, for = 250 GeV and tan/5 = 30. m^o is fixed to 250 GeV 
for the ohservahles Oh", Oh°, Oh+, and Oa°- For Ot, mA° = 100 GeV. 
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Figure 7: Predictions for the observables including the SUSY QCD correc- 
tions for MsusY = Mg = \fi\ and fi < 0. The solid contour lines follow 
the points in the {mA°,tsji (3) plane with constant value of O. The shaded 
area represents the region where the corrections are smaller than the theo- 
retical uncertainty. The long (short) dashed lines join the points where an 
experimental resolution of 50% (20%) is required to achieve a meaningful 
measurement. 
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Figure 9: Predictions for the observables with the SUSY QCD corrections 
reduced by 50% in order to simulate the effect of the SUSY EW contributions 
in the conservative scenario discussed in the text. The rest of the inputs and 
specifications are as in Fig. j^. 
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would approximately double the ones in Figs. |^ and |[ Also, the shaded area 
has doubled. The points with tan /? > 45 — 55 are testable at the one sigma 
level with a resolution of 50%; tan (3 > 20 — 25 with a resolution of 20%. In 
summary, even in this pessimistic case providing SUSY EW corrections of 
considerable size, a large region in the (m^,tan/?) plane remains testable. 

5 Conclusions 

In this paper we have looked for indirect signals of a heavy supersymmetric 
spectrum via its contributions to the radiative corrections to Higgs boson 
decays. For that purpose we have analyzed the dominant SUSY QCD cor- 
rections, coming from heavy squarks and gluinos, to the Higgs boson decays 
into quarks. We have studied in detail the nondecoupling contributions, es- 
pecially in the large tan (3 region where their size is large. In order to search 
for these SUSY signals, we have proposed a set of observables consisting of 
ratios of Higgs branching ratios into third-generation quarks {H°, h°, A° — >■ bb 
and tb) divided by the corresponding ones into third-generation lep- 

tons {H°,h°,A° T^T^ and z/r+). In addition, the observable for 

top quark decays given by the ratio of B(t H^b) divided by B{t W^b), 
complementary to the previous charged Higgs boson observable in the low 
mA° region, has been analyzed. These observables are optimal for this pur- 
pose since the SUSY QCD corrections appear just in the decays to quarks 
and, therefore, the decays into leptons can be used as control channels. 

We have carefully studied any sources of uncertainty that would modify 
the prediction of the observables, previous to the SUSY QCD corrections. In 
particular, the theoretical uncertainties coming from the QCD corrections, 
and from the errors in the values of the SM parameters involved in the 
determination of these observables, have been evaluated. In addition, it has 
been found that the SUSY QCD corrections will allow one to discriminate 
between the MSSM and a nonsupersymmetric 2HDMII, even in the case of 
a similar Higgs sector mass pattern. 

Our detailed study of this set of observables has revealed their high sen- 
sitivity to the SUSY QCD contributions and shown that they are sizable in 
most of the (m^o, tan/?) plane. The corrections to the different observables 
are strongly correlated. A global analysis of all of them, together with the 
experimental determination of m^o and tan /3, will be part of the search for 
supersymmetric signals at future colliders, especially if the SUSY spectrum 
turns out to be very heavy. The measurement of these observables should be 
the next step after a Higgs boson discovery at the LHC or the Tevatron. We 
have seen that, with a modest experimental resolution of 50%, the observ- 
ables would show evidence of the SUSY QCD corrections to the H°, and 
decay widths, if tan /3 > 20 — 25. A better resolution of 20% would show 
the corrections down to tan /? > 10 — 15. Our study of the SUSY EW effects 
shows that, even in the pessimistic case where the SUSY QCD corrections 
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get reduced by 50%, an experimental resolution of 20% would expose the 
SUSY corrections down to tan (3 > 20 — 25. Observation of the corrections 
in the h° case would require, in general, quite large values of tan f3. We hope 
our study will encourage our colleagues of CDF, DO, ATLAS, and CMS to 
investigate the actual performance of their detectors. 
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